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SUMMARY 


Three-Dimensional linear secondary instability theory is extended for 
compressible boundary layers on a flat plate in the presence of finite amplitude 
Tollmien-Schlichting (T-S) wave* The focus is on principal parametric resonance 
responsible for strong growth of subharmonics in low disturbance environment. 


I. INTRODUCTION 


Theoretically or experimentally, the compressible stability theory lacks 
a firm connection with boundary-layer transition. "There is little doubt 
that transition is preceded by linear instability in many instances, but the way 
these individual unstable waves. act, alone or in combination , to trigger the 
transition process is not known*' [1]. The nonlinear theories and the secondary 
instability that are much a prominent feature of incompressible stability theory 
do not exist for compressible boundary layer. 

Recently, Erlebacher and Hussaini, [2,3] by using direct simulation of 
parallel boundary layer, generated numerically a high Mach number vortical 
structure (peak-valley) similar to that observed and computed for incompressible 
flows . 

Stimulated by this work, we formulated the secondary three-dimensional 
instability problem for compressible boundary layers. In this paper, we 
investigate theoretically the effect of finite amplitude two-dimensional (2D) 
wave on the growth of three-dimensional (3D) perturbations in compressible 
boundary layers. Hence, this paper covers only a range of Mach numbers up to 
the transonic, where the critical primary disturbance is 2D. For supersonic 
boundary layers, the critical primary disturbance is 3D, and it is investigated 
by the author elsewhere. An analysis similar to that of Herbert [4], 
Herbert et al. [5] and Nayfeh [6] is followed here but spatial stability is 
considered for both primary and secondary instabilities. 

The primary instability leads to the growth of T-S waves and a streamwise 
almost periodic modulation of the flow. We study the linear stability of this 
flow with respect to spanwise periodic 3D disturbance. Floquet theory gives as 
a solution to the stability equations, all various types of resonance. We 
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consider the case of principle parametric resonance responsible for strong 
growth of subharmonics in low disturbance environment. 


II. PRIMARY INSTABILITY AND THE BASIC FLOW 


We consider the primary instability of the 2D compressible boundary layer on 
a flat plate, with respect to 2D quasi-parallel spatially growing T-S 
disturbances. The critical disturbance is 2D up to the transonic range of Mach 
numbers [7]. The flow field is governed by the 2D compressible Navier-Stokes 
equations , the energy equation and the state equation. Dimensionless quantities 
are introdu ced by using the reference velocity U r and th e reference length 
L = yUx/ Uf, so that Reynolds number is given by R = VU r x/i/ , where x measures 
the distance from the leading edge of the plate, and is the fluid kinematic 
viscosity. The thermodynamic and transport properties of air (treated as 
perfect gas) are made dimensionless using their corresponding freestream values. 

At sufficiently large distance from the leading edge, primary instability 
of the compressible laminar flow occurs with respect to T-S disturbances. These 
disturbances take the traveling wave form, 


(1) (u p ,v p ,p p ,t p ,p p ,p p ) = A o [ u(y),v(y),p(y),t(y),p(y),p(y) ] 

exp [i( j°<. dx -iot)] + c.c. 

Where the y-axis is defined normal to the flat plate, u and v are the 
disturbance velocities, p is the disturbance pressure, t is the disturbance 
temperature, p is the disturbance viscosity, and p is the disturbance density. 

For the spatial stability analysis o< is a complex wavenumber given by 
c* =c< ( . + ic*; and 6} is a real disturbance frequency, and c.c. denotes complex 
conjugate terms. The eigensolutions u,v,p and t are governed by a six-order 
system of equations that is given in reference 9. The density disturbance is 
related to the temperature and pressure disturbances through the state equation, 


(2) Y’Mp = pT + p t 


while the viscosity disturbance is related to the temperature disturbance by 

(3) p=(dp/dT)t = p t 

where p and p are the mean-flow density and viscosity respectively ,/ is the 
ratio of specific heats and M is the freestream Mach number. The six-order 
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system is numerically integrated as initial value problem using a combination of 
shooting [8] and Newton-Raphson iteration technique that employs a Gram-Schmidt 
orthonormalization procedure. The solution of this eigenvalue problem is 
calculated by the author elsewhere [7,9]. 


The linear stability theory of primary instability provides ©< for a given to 
and R. Then the integration of the growth rate -©<,• gives the amplification 
factor , 


(4) 


In ( A / A q 


f R 

-2 J c<j dR 
R op 


where A-, is an arbitrary initial amplitude of the primary instability at Rop (R 
where the onset of the primary wave) . The eigensolutions may be normalized 
such that A measures directly the maximum r.m. s. value of the streamwise 
disturbance, that is 


2 

(5) max | u(y ) J = 1/2 

0$ y < oo 


Since the primary instability of boundary -layer flows is induced by viscosity, 
the growth rates and amplification factors here are typically very small. 


Now the basic state under consideration is composed of the sum of the 2D 
compressible mean flow and a finite amplitude A of the primary T-S wave, that is 

(6) ( u b , v b ,P b ,t b ) = (U,0,P,T) + A [ u(y),v(y),p(y),t(y) ] exp(i0) 

Where 

A = A 0 exp(-J o<j dx), assumed constant, and 


9 - J c< r dx - W t 


III. SECONDARY INSTABILITY 


We consider the 3D quasi-parallel spatial subharmonic instability of the 
basic state given by (6). The finite amplitude primary wave acts as a 
parametric exitation on the secondary instability. Following the analysis of 
Herbert [4,5] , Nayfeh [6], and El-Hady [10] we apply Floquet theory and express 
the secondary wave using the normal mode concept. 
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(7) 


(u 


s > v s 


>P S »>-s 


»t s ) 


= exp(J Y dx)[u(y),v(y),p(y),t(y)] exp(^ i0) cos pz 


+C.C. 


A 

w s = exp(J Y dx)w(y) exp(^ i0) sin + c.c. 

A 

where p is a spanwise real wavenumber, and Y = Vp + i Yf is a characteristic 
exponent. The spatial growth rate of the secondary wave is given by Y^> while 
Yj can be interpreted as a shift in the streamwise wavenumber. In our 
calculations, we consider only the case of Y{ =0, that is the secondary wave is 
perfectly synchronized with the basic state. 

The secondary wave (7) is superposed on the basic state (6) and the result 
is substituted into the dimensionless Navier-Stokes equations. The meam flow 
plus the 2D T-S quantities are substracted, and the resulting equations are 
linearized in the secondary disturbances. Then one obtains an eigenvalue 
problem that can be written as 


-l -l 2 -l 

(8) I ( gZ, -T DTZj+ DZj + p ) + G ( Y M Z^ - T Z^ ) 

+ A (aj) =0 

- -I 2 2 -l 

(9) { G - p R [ (m+2) g - p ] } Z t - R D( p Z 2 ) 

- -I - -I 

+( T DU - Dp R g ) Z s - (m+1) p R g (DZ S + p Z^ ) 

~\ A A 2 

+ g Zg - R D(p DU) Z-j - p R DU Z ? + A (a^) + A (b 2 ) = 0 


(10) [ G - p r“' ( g + p 2 ) ] Z 3 - (m+1) p R~‘( g Z 2 + p Z g ) 


-1 ~ -l 

-(m+2) R D(p DZ^ ) - m R Dp(gZ, + pz^) + DZ fi 

a -I 2 

- p R DU g Z^ + A (a^) + A (b^) = 0 


I 


694 



(11) (nrt-1) p R* 1 J5 ( g Zj + DZ 3 ) + Dp R* l p Z g 

_ -I 2 - 2 -l 

[ G-fiR g + (m+2) p R p ] Z^ - R D( p Z g ) 

2 

- P Z 6 + A (a 4 ) + A (b^) = 0 


( 12 ) 



- -1 -» * 

p R DU ( Z-, + g Z^ )+ T DT Zj - M T G Zg 

-1-1-22 - a -i 2 A 

+{ G - R Pr [p (g - |3 ) + D(p DT)] - M R (DU) p } Z-j 

-l -l - - 2 

- R Pr ( 2 Z g Dp - p DZ ? ) + A (a fi ) + A (b g ) = 0 


where 


A A 

Zj = u , Z^ = Du 

A A 

Z g = Dw , Zg = p , 
g = ( + - i c*,- ) 

G - t" 1 ( gU - 2 i^) 


A 

M 


<Y- 


1 ) M 


A 




> 



A 

Dt 


also Pr is Prandtl number, D*d/dy, and m=2(e-l)/3 is the ratio of 
coefficients of viscosities, where e=0 corresponds to the Stokes hypothesis. 


The boundary conditions are 


(13) 


Z, - Z 3 = Z A = Z-j =0 at y- 0 


Zj , Z 3 , Zq , Z-j ■* 0 as y ® 


the 
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Equations (8)— (12) govern the secondary 3D instability of compressible 2D 
flows. They represent the mass, x-momentum, y-momentum, z-momentum, and energy 
equations respectively. The coefficients of A and A2 (a; , b( , i=l-5 ) are given 
in Appendix A. The coefficients of A in the mass equation and those of in 
all other equations are new in this system of equations in contrast to the 
secondary instability equations for incompressible flows. These coefficients 
are present mainly due to the density disturbance of both primary and secondary 
instabilities. In the incompressible limit equations (8)-(12) reduce exactly to 
that given by El-Hady [10]. When A=0,the system of equations (8)— (13) govern a 
primary subharmonic 3D wave. 

The system of equations (8)— (12) can be written as eight complex equations 
in the form, 


DZj - Z 2 = 0 

DZ 2 + c Q DZ Z + Cj = 0 

DZ S + c 2 DZ 3 + c 3 = 0 

DZ - Z- - 0 

(U) M 

DZ & + Cq DZg + Cg = 0 

DZ g + c 6 DZ 6 + C 1 = 0 

DZ-j - Z g = 0 

DZ ? + c y DZy + Cg = 0 

where the c's are quadratures in the primary and secondary disturbance 
quantities, and the overbar indicates a complex conjugate. Equations (14) are 
numerically integrated as initial value problem from y - y e (edge of the 
boundary layer) to the wall. The eigenvalue search used a Newton-Raphson 
iteration technique to satisfy the last boundary condition at the wall. A well 
tested code SUPPORT [8] is used which is coupled with an orthonormalization test 
based on the modified Gram-Schmidt procedure to overcome the stiffness of the 
integrated system of equations. 

The linear stability theory of the secondary instability provides for a 
given /J and R. Then the integration of the growth rate Y r gives the 
amplification factor, 

f R 

(15) In ( S / S Q ) - 2 ) V r dR 

Ros 
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where s D is an arbitrary initial amplitude of the secondary instability at R 0 s (R 
where the onset of the secondary wave) . The secondary subharmonic instability 
is believed to originate from a strong mechanism of combined tilting and 
stretching of the vortices such as the case of incompressible flows [10], thus 
the growth rates and amplification factors are expected to be large as they 
occur on a convective length scale. 


IV. RESULTS AND DISCUSSION 


For the incompressible limit (M = 0), our compressible secondary 
instability code gives results that are in full agreement with those obtained by 
Herbert et al.[5] and by El-Hady [10]. All results reported here are for the 
nondimensional frequency F * U> / R = 60E-6 , that remains fixed as a wave of 
fixed physical frequency travels downstream. 

At M * 0, a primary 2D instability grows between R op « 554 and R^ p = 1052 
(first and second neutral points) reaching a maximum amplification lactor of 
A/A0 * 41.68. As Mach number increases, the growth rates of the primary 2D 
waves decrease as shown in Fig.(l). Also the first and second neutral points, as 
well as the streamwise location where the maximum growth rate occurs^ shift to 
the left, they occur earlier upstream as Mach number increases. Fig. (2) shows a 
reduction in the amplification factor of the primary 2D instability as Mach 
number increases. 

Almost in the same region where the primary 2D waves are growing, a broad 
band of spanwise wavenumbers of primary 3D subharmonic waves are subject to 
amplification. Fig. (3) shows the growth rate curves of these primary 3D 
subharmonics (F * 30E-6) for different Mach numbers. These curves possess the 
same features of the primary 2D waves of having lower growth rates that shift to 
lower R as Mach number increases. Fig. (4) shows the amplification factors of 
these primary subharmonics. Both Figs. (3) and (4) suggest that the growth rates 
and amplification factors of the primary subharmonics are typically so small to 
bear any resemblance to experimentally observed transition. 

However, the growth rates shown in Fig. (5) as function of the spanwise 
wavenumber B -10E+3 p/R for secondary 3D subharmonics, are much larger than 
those for primary 2D waves or primary 3D subharmonics. This strong growth is 
due to the parametrical exitation by the finite amplitude primary 2D wave. 
Fig. (5) is calculated at R-850 and a primary 2D amplitude A *.01, for different 
Mach numbers. For high spanwise wavenumbers (B^.225), compressibility appears 
to have a stabilizing effect on the secondary subharmonic instability. In other 
words, the secondary subharmonic instability is largest at M * 0, and decreases 
as Mach number increases. But for low spanwise wavenumbers (B < .225) the effect 
of compressibility depends on the value of the spanwise wavenumber. 

At fixed F and R, results for different Mach numbers show a destabilizing 
effect ( higher growth rates for the secondary subharmonic ) as the primary 2D 
amplitude A increases. At R-1050, Fig. (6) illustrates the destabilizing effect 
of increasing A for M -0.8 compared with that for M -0. Fig. (6) reveals that 

the influence of compressibility on secondary subharmonics is not the only function 
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of the spanwise wavenumber, but also is a function of the primary 2D amplitude. 
While at A*. 01 compressibility is stabilizing for large values of B and 
destabilizing for small values of B, it has a destabilizing influence on the 
secondary subharmonics at A=.002 for all range of spanwise wavenumbers. The 
figure also shows that at very small amplitudes, considerable growth rates exist 
in a small band of wavenumbers, that extend to larger values as the amplitude 
increases • 

Fig. (7) shows the effect of Reynolds number R on the growth of the 
secondary subharmonics at M =0.8. As R increases , an increase in the growth 
rates exist at fixed F and A. 

To evaluate the overall effect of compressibility on the secondary 
subharmonics, we can combine the influence of increasing the amplitude A of the 
primary 2D wave and increasing R for various Mach numbers at fixed F. For 
comparison purposes, the amplification factor of the subharmonic is calculated 
using equation (15) from R os (onset of the secondary subharmonic instability) to 
any R downstream. For different Mach numbers, Fig. (8) shows the variation of 
the growth rates of the secondary subharmonics with R at a spanwise wavenumber 
B=.15. The initial primary 2D amplitude used for these calculations is A=.001. 
Fig. (9) shows the amplification factors decreasing from about 29 at M =0. to 12 
at M =0.8. 


Figs. (8) and (9) indicate that while the growth rates and amplification 
factors of the secondary subharmonics are decreased by increasing Mach 
number, the onset of the secondary instability is almost not affected. This is 
probably due to the combined effect of the upstream shift of the onset of the 
primary 2D instability, as well as the growing of the primary amplitude. 


V. CONCLUSION 


We formulated the secondary three-dimensional instability problem for 
compressible boundary layers. The effect of finite amplitude two-dimensional 
T-S wave on the growth of three-dimensional subharmonics is investigated for a 
range of Mach numbers from 0 to 1.2. Numerical results for F -60E-6 show that 
the local (at fixed R) effect of compressibility on the secondary subharmonics 
may be stabilizing or destabilizing depending upon their spanwise wavenumbers, 
as well as the finite amplitude of the primary 2D wave. However, the overall 
effect of Increasing Mach number is a reduction in the growth rates and 
amplification factors of the secondary subharmonics, almost with no change in 
the streamwise location where this instability sets in. 
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APPENDIX A 


a, = f(g Z, + p ) + D(f Zj) + g u F + D(v F) 
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b | * 0 


__l — — A 2 _ 

a., = [ G f + T gu-m 2 Rg(g + io<.)t + R p t]Zj 


— l — 1a — -\ -t _ a - 

+ [ T v - R D( [j t)] Z 2 + [ f DU + T Du - R g D( [i t)] Zj 


A _ ,»v _ _ A _ _ 

- R t DZ 2 - R ( mg + g ) t DZ 3 + R p ( m, g + im<* ) t Z^ 

A 2 ~ 

- R [m 2 °< (i g + <*) u + + m g) Dv + D u ] Z-j 


a A - 

- R ( Du + i*v ) D( p Z^) + [ i(<*U -«) + v DU ] F 


b., = [(g + i<x) f u + g f u ] Z ( + ( f v + f v ) Z 2 + ( ^ Du + f Du ) Z^ 


+ ( v Du + v Du ) F 


_ a 2 _ 2 

a s =[fG + T (gu + Dv)-R(g + i g* - /3 ) t ] 

“l — I a — a2— — 

+ [ T V - m z R D( p t) ] DZj - m 2 R t D Z 3 + i(<*U - ^ ) v F 


_\ -l A — . — A _ _ -\ A - 

+ [IT e* v - m R D( p t) g ] Z ( - R (m, g + i«) t Z^ - m R p D( p t) Z^ 

A — A _ Z 

- m t R p t Zg - { R [ (lm<X+ g) Du + (i g -<*) v + m z D v ] 


-l A — a — 

- R Dp ( ii o< u + m z Dv )} Z-j - R ( im e* u + Dv ] Z y 


bj • 1<H( f v - f v ) Z, + [( f “ + f u ) g + f Dv + f Dv ] Zj 


+ (fv + fv) DZg + [ 2v Dv + ix( vu-vu )] F 
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- -I _ a _2 

a 4*t fG + T 8 U-R (g + ic<g-m fi ) t ] Zq 

-I -I * - * - 

+ [T v-R D( (j t) ] Z g - R [ t DZ s -p t ( g + m g ) Z, 


-I A - 

- m k p t DZ 3 - m p ( i <x u + Dv ) ] + R p D( p t) Z 3 


b 4 =(fu + fu)gZ^+(fv + fv)Z 5 


— - 1 _ A A _1 A _ 2 — 2 2 

a & = {f G + T g u - 2M R DU (Du + i <X v) - Pr [R t ( g + 2i<Xg - <* - p ) 


+ R - ' D( p Dt ) ] } i 7+ [ i (ex U -U> ) t + v DT ] F - M g u Z g 


- M v DZ g + [ T* 1 v - Pr”' ( 2 R Dt + t R Dp ) ] Z ? - R Pr ' t DZg, 


A A -I - 

+ [ i &< ( T t-Mp)-2HpR g ( m 2 u - m Dv ) ] Zj 


- 2 M r"* [ p (iexv + Du)+ p t DU ] ( Z 2 + g Z^ ) 


-I a - 

+ ( T Dt + f DT - M Dp ) Z 3 


a -l - - 

- 2M R p [ ( im <X u + m 2 Dv ) DZ^ + mp(i<*u + Dv)Z^] 


^ A 2 — • — . mm 2 ^ 

bg={ (fu + fu)g-2MR[©< vv + i«(vDu-vDu) + Oi u u 

2 _ _ _ 

+ m©<(uv+uv) + Dv Dv + Du Du ] 


+ [io<(ut-ut)+vDt + vDt]F + ( f v + f v ) Zy 


+ [ix(tf-tf)-2MR i g ( t u - t u ) ] Z ( 


- 2 M R ( t Du + t Du + io<( t v - t v ) (Z-^+gZ^) 


— — A A _ _ 

+ ( f Dt + f Dt ) Z^ - 2M R m 2 ( t Dv + t Dv ) DZ^ 
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Fig* 1 Growth rates of the primary 2D wave Fig. 2 Amplification factors of the primary 
(F=60E-6) at different Mach numbers* 2D wave (F=60E-6) at different Mach 

numbers . 
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Fig. 3 Growth rates of the primary 3D Fig. 4 Amplification factors of the primary 
subharmonic wave (F=30E-6 and B=.12) at 3D subharmonic wave (F=30E-6 and B=.12) at 
different Mach numbers. different Mach numbers. 
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Fig. 5 Growth rates of the secondary 3D subharmonic at 
R=850 and A=.01 as function of the spanwise wavenumber, 
for different Mach numbers. 
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Fig. 6 Effect of the primary amplitude A on the growth 
rates of the secondary 3D subharmonic at R=1050, for 
Mach numbers 0 and 0.8. 
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Fig. 7 Effect of Reynolds number R on the growth rates 
of the secondary 3D subharmonics at M =0.8, for 
different primary amplitude A. 
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Fig. 8 Variation of the growth rates with R Fig. 9 Variation of the amplification 
of a secondary 3D subharmonic at B=.15 for factors with R of a secondary 3D 
different Mach numbers. subharmonic at B=.15 for different Mach 

numbers. 
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